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ABSTRACT 

We carry out next-to-leading order (NLO) parton model calculations for the 
standard hard "QCD" processes in the massless Schwinger model. The asymp- 
totic expansion of the exact result for the deep inelastic cross section is used to 
infer the NLO distribution function. These distribution functions are then used 
to calculate the NLO Drell-Yan parton model cross section and it is found to 
agree with the corresponding term in the expansion of the exact result for the 
Drell-Yan cross section. Finally, by using the bosonization formula and the exact 
solutions we study the interference between different partonic processes. 
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1 Introduction 



Although the parton model |1|, ^ lacks a solid theoretical foundation, it is an essential tool 
in QCD applications to strong interactions and so far it has been reasonably successful. If 
we have a theory that exhibits many properties of 4-dimensional quantum chromodynamics, 
QCD4, and this theory is exactly solvable, then it is important that we study how well parton 
model works in this theory. The massless Schwinger model is one such theory in that it is 
an exactly solvable, interacting quantum field theory P, ^ that is both asymptotically free 
and exhibits confinement p, 

In reference 0, the exact cross sections for lepton-antilepton annihilation, deep inelastic 
scattering and Drell-Yan processes in the massless Schwinger model coupled to a scalar 
current are calculated in terms of the functions R±{q^)- In the full asymptotic expansions 
of R±{q^) were described and the terms up to and including order {g'^/q^Y were explicitly 
computed where g is the strong coupling constant and is the squared momentum transfer. 
The leading terms in these exact cross sections were then shown to equal the leading order 
parton model results in the Bjorken scaling region. This paper extends this analysis of the 
massless Schwinger model to the next-to- leading order (NLO). 

Of course, the lack of transverse momenta for the Schwinger model partons and the corre- 
sponding absence of jets limits the extent of the lessons that we adduce for QCD4. Another 
fundamental difference between the Schwinger model and QCD4 is that the Schwinger model 
coupling constant g has the dimensions of mass. However, we see this difference as an op- 
portunity to consider several potential problems for the parton model when it is is pushed 
to next order in g"^ /q^ including: (1) the mass of the hadron which is of order g in this model 
can no longer be neglected in the kinematics; (2) hadronization or bound state effects are as- 
sumed to be suppressed by inverse powers of q^ and so may arise at order g"^ /q^; (3) quantum 
interference of the hard partonic processes are normally assumed to be power suppressed and 
so may arise at order g"^ /q^; and finally (4) the analogue of higher twist terms. 

To our surprise, despite these complications, we find that the parton model NLO Drell-Yan 
cross section agrees with the exact Drell-Yan cross section at NLO. We also argue that at 
NLO the Schwinger model photon, the analogue in our model of the gluon, has a parton 
distribution, which may be calculated by evaluating processes at order g^/q^. As another 
application, we use the bosonization formula to isolate the contributions coming from the 
interference of the underlying partonic processes. Then from the expansion of the exact 
results, we find that the interference is suppressed by order {g"^ /q^Y- 

2 NLO Parton Distribution Functions From Deep In- 
elastic Scattering 

The processes that we calculate occur in the same model used in [0, ^ ^ ^. For completeness 
we briefly describe the model here. We extend the massless Schwinger model by including 
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a massless fermion / (our "lepton" ) that is not QED2-charged but interacts with the QEDg- 
charged fermions, through a Yukawa couphng to a scalar photon (p. The full Lagrangian 
is §,§ 

L = FV4 + ^{i^ + g m + n + l/20n0 + e(V>V' + /7)0. (1) 

So e is the analogue of the electromagnetic coupling, and g is the analogue of the QCD4 
strong coupling constant so from now on, we will call the photon of QED2 as "gluon". All 
of our calculations are done at lowest order in e. Then, using the dual realization of the 
Schwinger model in terms of a free scalar of mass( which we will call "hadron") rrih = g/\^ 
and the bosonization formula 

tjjtjj = c: cos{2y/TTh) :, (2) 
we may compute the process cross sections exactly in g. 



First consider the DIS process 

f(k,) + hiP)^fik2)+X, 

where h is the target particle and the final states X are summed over and where the mo- 
mentum assignments are given in the parantheses. The cross section at leading order in e^, 
but all orders in g"^, is 

da = -^-^R,{{q + Pr)dk'. (3) 

where g = fci — ^2 is the momentum transfer, and P is the hadron momentum, R^lq"^) is 
defined by 

R,iq') = {R+{q')-R-{q'))/^ (4) 

and 

R±{(f) = (^j 'i^a;exp(zgx) exp(±47rAm(x)), (5) 
In 1^ we showed that for 7^ 0, 

R^{q^) = l + + (^2/^2)3 J_(i2 + 4 ln[V/g'])+0((//g')'), (6) 



and 



R-{q') = {gVq'r^{6.96 + 3.79H7rgyq']) + oiigyq'r). (7) 



Thus the cross section in DIS has no o{g'^ /q^) correction to the leading order result. We will 
use this and the parton cross sections presented next to evaluate the NLO parton distribution 
functions. 



*The prefactor c is a normal ordering dependent constant and it equals c = (2^Y /2 '^hen the normal 
ordering mass equals m/j. 
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The parton model calculation of DIS is done perturbatively in g where the final states 
that are summed over consist of quanta of ip and of A^^. The vector field A^^ has asymptotic 
propagating states because we have chosen to regulate the IR collinear divergences (the same 
fermion mass singularities that occur in perturbative QCD4 calculations) by temporarily 
giving the A^^ field a mass, rrig. At leading order in g the only process that contributes is 

/(fci) + ^(Pi)^/(A;2)+^(P2), 
and the cross section is given in to be 

where z = —q^/(2pi ■ q). 

As in QCD4, the NLO g^ corrections come from the interference of the one loop corrections 
with the leading order process, dal, and from gluon bremstrahlung, dal- We choose to 
regularize by giving the "gluon" a mass uig. The calculational steps that go into evaluating 
da]j are almost the same as in those outlined in |P. The result, written in terms of /3 = m^/g^, 
is found to be 



The calculation of the bremstrahlung^ cross section is a little more involved and warrants a 
brief description. The process is 

/(fci) + ^(Pi) ^ fih) + V^(P2) + A^{p,), 
and the corresponding cross section is directly calculated to be 

,.1 dk2g'^ 1 (-8g^) dp2dp3 ^ ^ nn^ 

d^b = Tl~^ — 2^7 ^^7 ■ P^~^~^^ ^1 +Pi-k2-P2- P3 • 10 

4g4 Ep q^ 2tt {pi - pzY{p2 - PzY E2 E^ 

It is an exercise in kinematics to do the integrals over p2, Ps and rewrite this cross section 
in terms of the two independent variables and z. The result is found to be 



I dk2 g"^ I 1 



'111 



For this radiation correction it is important to note that the condition 

(pi + qf > ml, (12) 



^Throughout this paper, we denote the cross-sections of the partonic subprocess, the parton model cal- 
culation, and the exact calculation by a, a, and a respectively. 

^2-diniension gauge bosons have no physical degree of freedom, however, after regulariztion, the Schwinger 
model phonton obtains a longetudinal mode. This makes the model much more like QCD. This method has 
proven successful in the calculation of lepton-antilepton annihilation process. (Q) 
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implies 

Now we are in a position to fold the partonic cross sections into a hadronic result. First 
notice that the mass of the hadron is of order g, so we cannot neglect this mass in the 
next-to-leading order calculation. We will follow the well-established method to deal with 
the hadron mass^]. We define, without loss of generality, the hadron momentum fraction 



y carried by the parton as 

pi+ = yP+, (14) 

where 

pI = o , q+ = qo + qi , q- = qo-qi- 

It follows that z = ^/y, where ^ is the Nachtmann scaling variable 

2r 

e = , (15) 

1 + Y^l - 4:X^ml/q'^ 

Also the familiar Bjorken scaling variable x is defined to be 

The hadronic cross section in the parton model is then given by multiplying the underlying 
partonic cross sections by the parton distribution functions, both given as a functions of ^, 
y, and g^, and integrating over y over its allowed limits. Thus 

da = J2 f'dyMy,q')da{y,q'). (17) 



The partonic cross sections in equations (§),(0), and (0) must be summed to give da to 
the desired order. Further charge conjugation invariance of the underlying theory implies 
that the fermion and antifermion distributions are equal. Hence the sum over i in (|T7p is 
replaced by a factor of 2. Finally, the leading order result / = 1 is all that is needed in 
evaluating the integral over y with the NLO cross sections dal and dal- The only subtlety 
in the integration over y occurs in the integration limits for the bremstrahlung cross section. 
There equation (|^) implies that 

(l-/3)e<2/<l. (18) 

Setting the partonic result equal to the exact result and keeping terms up to order g'^/q^ 
gives 

/(2/,?') = l-4^1n[^]. (19) 
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This is the full NLO distribution function for Schwinger model partons in the DIS process. 
Recall that in QCD4 the NLO distribution functions may be characterized in terms of the 
Altarelli-Parisi equations [|n|. These differential equations describe the running of the NLO 



distribution functions with q due to the collinear divergence of the radiative gluon correc- 
tions. The initial conditions for the distribution functions must come from a comparison with 
experimental data at some fixed q^. In the calculations of this section, we have effectively 
solved for both the running and the initial conditions by directly relating the cross sections. 

Observe that the momentum sum rule is not satisfied by these NLO distributions. Indeed, 

2 /' dyyf{y, q^) = I + ^ x (-00). (20) 
Jo rrig 

From this we conclude that, in the massive QED2-gluon regularization scheme, the gluon 
has a distribution whose expansion in g starts at order g'^/q'^. To calculate this distribution 
from DIS would require going to order g^/q^ since the contributing cross section is itself of 
order g'^/q^. Instead, we will now use the NLO distribution functions inferred from DIS to 
evaluate the Drell-Yan cross section. 



3 Drell-Yan Process at Next-to-leading Order 

We next use the parton model to calculate the cross-section for the Drell-Yan process, 

h{P) + h{P')^ f{k)+f{k')+X. 

In the parton model it is assumed that the distribution functions /(x, q^) used in the calcu- 
lation of DIS are the same as the ones /^^(x, q^) to be used in the calculation of Drell-Yan 
process [0. But, a prior, these distribution functions /^^(x, g^) with timelike momentum 
transfer, q^ > 0, and f{x,q'^) for the DIS process with spacelike momentum tranfer, q^ < 0, 
could be different. In [0 it is shown that they are both equal to 1. However, higher order 
radiation might spoil this equality. In fact, the exchange of soft gluons between the initial 
state hadrons is often cited as a possible complication in even applying the parton model 
to the Drell-Yan process. Nevertheless, we will verify in this section, that the NLO parton 
model cross section, with the initial state radiation included, equals the exact cross section 
at NLO. 

The cross section to leading order in but exact in g'^ is given in to be 

p'^ 1 

da=— rdkdk'RUq- P - P'f), (21) 

q'[{p.p')2-gy^2^l 

where 

Riq') = R+{q') + R-{q') 

= J d'^x exp{iqx) cosh{ATcAm{x)). (22) 
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Then using the results of the asymptotic expansion in equations (^|^ we see that the exact 
Drell-Yan cross section is 

^^ = ?p:i^^^ + ^V^)- ^^^^ 



To evaluate the NLO parton model, Drell-Yan cross section we must evaluate the underlying 
partonic cross sections and then fold them together with the parton distributions in the 
incident hadrons. At leading order the parton process is 

and the cross section is evaluated in to be 

da%Y = ^dkdk'S\p + p' -k-k'). (24) 
Zq 

To go to NLO we need to evaluate the radiative corrections to the leading order cross section. 
As with DIS, there are two corrections. One from the interference of the one-loop amplitudes 
with the leading order process, daj^yv- And the other from the emmission of a final state 
gluon from one of the incoming partons, daj^yb- 

The calculation of the virtual cross section follows the same lines as the calculations of the 
virtual corrections to DIS - given in the preceding section - and of the virtual corrections to 
/ — / annihilation described in M . The result is that 



dcr},y^, = da%y^-^^ln\(3\. (25) 

(J Tip 



The other order g'^/q'^ correction comes from the process 

m+^{p')^f{k) + f{k') + A^{r), 
and the corresponding cross section is immediately calculated to be 

d^'DY, = —,dkdk'5\k + k' + r-p- (26) 

These partonic cross sections are then folded together with the NLO distribution functions 
to get the NLO parton model cross-section as follows 

da = YlJ d,ydy'fi(y)fj{y')da,j{y,y')dkdk' 

= 2 dydy'f{y)f{y')[daly + da],y^^ + j dr-^]. (27) 

In the second line of equality, we have used charge conjugation invariance to replace the 
sum over parton species with a factor of 2. Also, the outgoing gluon momentum must be 
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integrated over in the bremstrahlung contribution. We will evaluate the contributions from 
the three cross sections separately. 



For the first two terms: 



2 / dydy'f{y)f{y') [da^y + da\ 



DY,v 



—dkdk' / dydy'S'^iq -p-p'){l , a 
Jo I np 



Inl^ + ln'-^' 



y 



y' 



+ 4^1n/?| (28) 
q^ Tip I 



where 



p+ = yP+ , p'_ = y'P'_ . 
The integration over y and y' can be easily done using the following identity: 

6'{q) = 25(g+)5(g_) 

and we get: 

2 dydy'f{y)f{y') {da%y + da^y^^) 



1 

qi pp, 



dkdk' 



1 , (2PP' -2Pq-2P'q + q^^ 
1 ^ — ^ In 



q^ TiP \ 



(29) 



(30) 



Finally, consider the contribution due to the bremstrahlung cross section, daj^yf^. We need 
only the leading order distribution function, and at this order, we can neglect the effect of 
the mass of the hadron, but now the integral is complicated by the extra integral over the 
gluon momemtum, r^. Using equation (|29|) , we can do the integrals over y and y' first, and 
after some algebras, get: 



^^^^^ ' dkdk'C^ft 



2 / dydy'f{y)f{y') dr ^ , ^, ^ 

/o J dr q^ PP' q^ nfj J E., 



(31) 



We are left with the integral over the momentum r. The limits of the integration are 
determined by the conditions that y < 1 and y' < 1 which imply that 

2(P:-g_) - - 2(P+-g+) ^ ^ 

Therefore, the contribution to the parton model Drell-Yan cross section due to gluon emission 
is given by 



2 / dydy'f{y)f{y') I dr 
Jo J dr 



a 



1 

DY,b 



1 ^,^,,(7^ 1 , f2PP'-2Pq-2P'q + q' 



-dkdk — — -In , 
q4 pp, q2 y 



q'(3 



(33) 



Comparing equations (|30D and (|33D, we find that the infinity in the parton distribution 
function cancels with that in Drell-Yan parton cross section, and the net NLO correction to 
the Drell-Yan cross-section in the parton model calculation vanishes: 

.4 1 



daoY = daoY = -^^^dkdk' (34) 

That is, the NLO parton model calculation of the Drell-Yan cross section, using the NLO 
distribution functions inferred from the DIS process, equals the exact Drell-Yan cross section 
evaluated to next-to-leading order. 
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4 Discussion and Conclusion 



Simple dimension counting in QCD^ distinguishes tlie perturbative corrections according to 
their twist. Typically, one focuses on the perturbative large logarithmic corrections and 
neglects the higher twist corrections that are suppressed by powers of l/q^. In contrast, the 
solvable 1 + 1-dimensional massless Schwinger model coupling constant g has the dimensions 
of mass so that all of the corrections are suppressed by powers of jq^ ■ Thus the perturbative 
corrections can not be cleanly separated from corrections due to hadron mass, interference 
and hadronization, as well as potential higher twist terms. We have exploited this difference 
to study these effects in the parton model. 

It turns out the ^ scaling handles the hadron mass correctly in this model. However, it is 
kind of a surprise that we get a Drell-Yan cross section which agrees precisely with the exact 
solution without introducing multi-parton densities. There is still much to understand in 
this model. 

The calculations in this paper establish the universality of the distribution functions in 
the Schwinger model at next-to-leading order. Using the asymptotic expansion of the ex- 
act results we can test whether this process-independence persists to even higher orders. 
Other assumptions of the model can also be tested. For example, in |§] we argued that the 
(gf^/g^)^ ln[(7^/g^] term in the expansion of the exact annihilation cross section was an effect 
of hadronization that can not be calculated from perturbation theory (see equation (|^). In 
QCD4 it is tacitly assumed that these hadronization effects are supressed by inverse powers 
of the large squared-momentum transfer. Another fundamental assumption of the parton 
model is that the interference of different partonic processes is suppressed by inverse powers 
of large squared-momentum transfer. For example, this assumption is built into the starting 
point of those proofs of various factorization theorems that rely on Landau-Cutovsky cut 
diagrams [|13|. At least in one case, we can test this assumption in the Schwinger model. 

First, separate the density in equation (||) that couples to the scalar photon into two parts: 



Let us use a and h to represent the annihilation operators for the quarks and antiquarks, 
respectively. Then 



Next, without loss of generality, assume the hadron is moving right. Then the contribution 
of the first term of equation ( ^5]) to DIS process corresponds to a quark coming out of the 
hadron and the second term corresponds to an antiquark coming out the hadron. In parton 
model, the cross sections for these two process are summed over. In the exact solution, 
however, we can calculate the interference of these two terms using 



(35) 



V'iV'fi ~ al^fi + ^L^L H 

^i?^L ~ flif flL + &L&iJ H • 



(36) 



(37) 
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Thus, it turns out the inference term is proportional to R-{q^) and from equation (0) the 
leading order of this interference occurs at o{g^/q^). This term is indeed is very small but is 
nonvanishing. A similar analysis applies to Drell-Yan process. 
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